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@ Fluid flow measurements are made in situ using a repeat formation tester witii a modified probe 
aperture, or on rock samples using a mini-permeameter with a modified probe aperture. The modified 
probe aperture has an elongate cross-section, such as elliptic or rectangular. A first flow measurement 
is made witti ttie longer dimension of the probe aperture in a first orientation {e,g„ horizontal or vertical) 
with respect to the fomiation bedding planes. A second flow measurement is made with the probe 
aperture orttiogonal to the first orientation, or witti a probe aperture of non-elongate (e.p.. circular) 
cross-section. Simultaneous equations relating values of known and measured quantities are solved to 
obtain estimates of local horizontal and/or vertical formation penneabOity. 
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BACKGROUND OF THE INVENTION 

1. Field of the Invention 

5 The invention concerns methods for estimating the horizontal and/or vertical components of permeability 
of an anisotropic earth formation. 

2. The Prior Art 

10 The permeability of an earth formation containing valuable resources such as liquid or gaseous hydrocar- 
bons is a parameter of major significance to their economic production. These resources can be located by 
borehole logging to measure such parameters as the resistivity and porosity of the formation in the vicinity of 
a borehole traversing the formation. Such measurements enable porous zones to be identified and their water 
saturation (percentage of pore space occupied by water) to be estimated. A value of water saturation signifi- 

15 cantly less than one is taken as being indicative of the presence of hydrocarbons, and may also be used to 
estimate their quantity. However, this information alone is not necessarily adequate for a decision on whether 
the hydrocarbons are economically producible. The pore spaces containing the hydrocarbons may be isolated 
or only slightly interconnected, in which case the hydrocarbons will be unable to flow through the formation to 
the borehole. The ease with which fluids can flow through the formation, the permeability, should preferably 

20 exceed some threshold value to assure the economic feasibility of turning the borehole into a producing well. 
This threshold value may vary depending on such characteristics as the viscosity of the fluid. For example a 
highly viscous oil will not flow easily in low permeability conditions and if water injection Is to be used to promote 
production there may be a risk of premature water breakthrough at the producing well. 

The permeability of a formation is not necessarily isotropic. In particular, the permeability of sedimentary 

25 rock in a generally horizontal direction (parallel to bedding planes of the rock) may be different from, and typ- 
ically greater than, the value for low in a generally vertical direction. This frequently arises from alternating hor- 
izontal layers consisting of large and small size formation particles such as different sized sand grains or day. 
Where the permeability is strongly anisotropic, determining the existence and degree of the anisotropy Is Inrv 
portant to economic production of hydrocarbons. 

30 Techniques for estimating formation permeability are known. One technique involves measurements made 
with a repeat formation testing tool of the type described in U.S. Patents No. 3,780,575 to Urbanosky and 
3,952,588 to Whitten, such as the Schiumberger RFT*^ tool. A tool of this type provides the capability for re- 
peatedly taking two successive "pretest" samples at different flowrates from a formation via a single probe in- 
serted into a borehole wall and having an aperture of circular cross-sectk)n. The fluid pressure is monitored 

35 and recorded throughout the sample extraction period and for a period of time thereafter. Analysis of the pres- 
sure variations with time during the sample extractions ("draw-down") and the subsequent return to initial con- 
ditions ("build-up") enables a value for an effective formation permeability to be derived for each of the draw- 
down and build-up phases of operation. 

Figure 1 illustrates schematically the principal elements of a tool employed in taking "pretest" samples. The 

40 tip 110 of a probe is inserted through mud cake 112 into the borehole wall. Mud cake 112 and a packer 114 
hydraulically seal the probe tip 110 with respect to the formation 116. The probe includes a filter 118 disposed 
in the probe aperture and a filter-cleaning piston 120. The pretest system comprises chambers 122 and 124 
and associated pistons 126 and 128. Pistons 126 and 128 are retracted in sequence each time the probe is 
set. Piston 1 26 is withdrawn first, drawing in formation fluid at a flow rate of, for example, 50 ccMiin. Then piston 

45 128 is withdrawn, causing a flow rate of, for example, 1 25 cc/min. Figure 1 shows the system in mid-sequence, 
with piston 126 withdrawn. A strain gauge sensor 132 measures pressure in line 134 continuously during the 
sequence. When the probe is retracted, the pistons 1 26 and 128 are moved to expel the fluid, and filter cleaning 
piston 120 pushes debris from the probe. 

The pressure measurement is recorded continuously in analog and/or digital form. Figure 2 shows a typical 

50 analog pressure recording during pretest. A pressure draw-down Api is recorded as piston 126 is withdrawn 
during a time period Ti, and a pressure draw-down Ap2 Is recorded as piston 128 is withdrawn during a the 
period T2. When pretest chambere 122 and 124 are full (at time i^, the pressure begins to build up over a time 
period At toward a final pressure, that of the formation. 

The permeability has been estimated by analyzing the pressure recording during either buildup or draw- 

55 down. As illustrated in Figure 3, the point 310 at which the probe tip 110 is applied to the wall of the borehole 
312 coincides with the center of the latter stage of the pressure disturbance during buildup. From the perspec- 
tive of a coordinate system whose axes have been suitably stretched by an amount dictated by the horizontal 
and vertical components of the permeabflity, the pressure disturbance appeare to be propagating spherically 
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outward from the probe tip 110. Thus the analysis yields a single "spherical' permeability value, consisting of 
a specific combination of both the horizontal and vertical components of the permeability. During drawdown, 
the pressure disturbance has only been analysed for the case of a homogeneous formation with isotropic per- 
meability. For the anisotropic case, the ad hoc assumption has been made that the isotropic permeability be 
5 replaced by the "spherical" permeability. Only in some cases could the analysis yield separate values for hor- 
izontal and vertical permeabilities, and then only with the Incorporation of data from other logging tools or from 
laboratory analysts of formation core samples. UntO recently, it had been assumed impossible to derive sepa- 
rate horizontal and vertical permeability values solely from the measurements provided by the single-probe 
type of tool. 

10 Another method of estimating formation permeability Is described in U.S. Patent No. 4,742.459 to l-asseter. 
Figure 4 shows in schematic form a borehole logging device 400 useful in practicing the method. In this ap- 
proach, formation pressure responses vs. time are measured at two observation probes (402 and 404) of cir- 
cular cross-section as a transient pressure disturbance is established in the formation 406 surrounding the 
borehole 408 by means of a "source" probe 410. The observation probes are spaced apart in the borehole, 

15 probe 404 (the "horizontal" probe) being displaced from source probe 410 in the lateral direction and probe 
402 (the "vertical" probe) being displaced from source probe 410 in the longitudinal direction. Hydraulic prop- 
erties of the surrounding formation, such as values of permeability and hydraulic anisotropy, are derived from 
the measured pressure responses. 

While the technique of this patent has advantages, the use of multiple spaced-apart probes has some in- 

20 herent drawbacks. Forexample, the MRTT^and MDT^tools commercialized by Schlumberger and employing 
principles of the Lasseter patent have the observation probes spaced some 70 cm apart along the borehole. 
The estimate of vertical permeability Is thus based on flow over a relatively large vertical distance. While this 
is sometimes appropriate, it is often preferable to obtain a more localized value of vertical permeability. If the 
longitudinally-spaced observation probes are set so that they straddle a hydraulic barrier In the formation (e.g., 

25 a formation layer of low permeability relative to the layers in which the probes are set), the values determined 
for vertical permeability and hydraulic anisotropy may differ significantly from the local characteristics of the 
formation layers above and below the barrier. Moreover, the technique of the Lasseter patent may require si- 
multaneous hydraulic seating of three probes, though it may be possible to make both horizontal and vertical 
measurements with only two probes. Accurate measurement may be prevented if one or more of the probes 

30 fails to seal properly, such as where the borehole surface is uneven. WhOe even a single-probe system can 
encounter seating problems, the need for simultaneous seating of multiple probes may increase the difficulty 
of obtaining the desired measurement. 

A method for determining the various components of the permeability of an anisotropic formation with a 
single probe is described in U.S. Patent No. 4.890.487 to Dussan V. etaL See also E.B. DUSSAN V. etal,, An 

35 Analysis of the Prassure Response of A Single-Probe formation Tester, SPE Paper No. 16801, presented at 
the 62nd Annual Technical Conference and Exhibition of the Society of Petroleum Engineers (1987). Pressure 
draw-down and build-up are measured as fluid samples are extracted from the formation at controlled flow rates 
with a logging tool having a single extraction probe of circular cross-section. This may be done with a system 
as shown in Figure 1, producing a pressure recording as shown in Figure 2. The measured build-up and draw- 

40 down data are analyzed to derive separate values for horizontal and vertical formation permeability. This is 
possible because they successfully analyze the pressure disturbance during draw-down for an anisotropic for- 
mation. This technique offers a localized determination of hydraulic anisotropy, and avoids the need to incor- 
porate data from other logging tools or core analysis. It has the disadvantage that it relies on measurement of 
pressure build-up, which demands an extremely fast-responding pressure transducer with a very high sensi- 

45 tivity. Pressure draw-down is a relatively robust measurement - pressure is measured before and after the 
pressure disturbance caused by fluid extraction. Pressure build-up is a more delicate measurement because 
the rate of pressure recovery must be measured accurately as the detected pressure asymptotically approach- 
es formation pressure (the pressure recovere at a rate of 1/t^). 

A further technique for determining permeability is performed in the laboratory using formatbn samples 

50 and a laboratory instrument known as a minl-permeameter. The instrument has an injection probe with a nozzle 
of circular cross-section which is pressed againstthe surface of a sample and appropriately sealed. Pressurized 
gas flows through the injection nozzle Into the rock sample as gas flow and Injection pressure are measured. 
Referring to the schematic view of Figure 5, the process maw be performed on a f iret face 510 having Its long- 
itudinal (z) axis perpendicular to the bedding planes of a formation sample 500 and on a second face 520 having 

55 its longitudinal (x or y) axis parallel to the bedding planes of the formation sample. The measured lows through 
the sample are used in determining permeability. See, for example, R. EIJPE et al., Geological NoteiMinl- 
Permeameters for Consolidated Rock and Unconsolidated Sand, THE AMERICAN ASSOCIATION OF PET- 
ROLEUM GEOLOGISTS BULLETIN. Vol. 55, No. 2, pp. 307-309 (1971); C. MCPHEE, PROPOSED MINI- 
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PERMEATER EVALUATION REPORT. Edinburgh Petroleum Equipment. Ltd., Edinburgh. Scotland (1987); and 
D. GOGGIN el a/.. A Theoretical and Experimental Analysis ofMinlpermeameter Response Including Gas Slip- 
page and High Velocity How Effects, IN SITU. 12(1&2). pp. 79-116 ( 1988). 

Determining horizontal and/or vertical permeabilities of a formation with a mini-permeameter has a number 

5 of important limitations. The mini-permeameter is a laboratory instrument, and cannot be used to make in situ 
measurements in a well bore. Thus, it can only be used to make the necessary measurements If formation core 
samples are available, which is not always the case. Moreover. It entails destruction of portions of the core 
sample, as a smaller sample having a smooth face parallel to and perpendicular to the bedding planes must 
be cut from the sample for testing. Also, the mini-permeameter measures the permeability of isotropic samples. 

10 In the case of an anisotropic sample; it only gives an effective value. Thus, it would only give an etfective vertical 
and effective horizontal permeability from the two faces 510 and 520, respectively. 

SUMMARY OF THE INVENTION 

15 It is an object of this invention to provide improved methods for determining horizontal and vertical per- 
meabilities of an earth formation. It is further an object of the present invention to proride methods which may 
be performed in situ or at the earth's surface. Another object of the invention Is to provide methods whteh avokJ 
limitations of the prior art methods described above. These and other objects are attained in accordance with 
exemplary embodiments of the invention described below. 

20 In a preferred embodiment, fluid low measurennenfs are made in situ using a repeat formation tester with 
a modified probe aperture, or a mini-permeameter with a modified probe aperture. The modified probe aper- 
ture has an elongate cross-section, such as elliptic or rectangular. A first flow measurement is made with the 
longer dimension of the probe aperture in a first orientation (e.g. , horizontal or vertical) with respect to the for- 
mation bedding planes. A second low measurement is made with the probe aperture orthogonal to the first ori- 

25 entation. or with a probe aperture of non-elongate (e.g.. circular) cross-sectton. Simultaneous equattons relat- 
ing values of known and measured quantities are solved to obtain estlnnafes of local horizontal and/or vertical 
formation permeability. 

BRIEF DESCRIPTION OF THE DRAWING 

30 

Preferred embodiments of the invention are described in more detaS below with reference to the aocom- 
panving drawing, in which: 

Figure 1 illustrates schematically the principal elements of a prior-art tool employed in taking "priest" for- 
mation fluid samples in a borehole; 
35 Figure 2 shows a typical analog pressure recording made during pretest sampling with a tool of the type 
shown in Figure 1; 

Figure 3 illustrates a prior-art model of a pressure disturbance in a formation; 

Figure 4 illustrates schematically a prior-art borehole logging device having a source probe and a spaced- 
apart pair of observation probes for formation testing; 
40 Figure 5 illustrates a formation sample used for mini-permeameter testing in accordance with the prior art; 
Figure 6 illustrates generally vertical fluid flow Into a horizontally-oriented, elongate probe aperture In ac- 
cordance with the invention; 

Figure 7 illustrates generally horizontal fluid flow into a vertically-oriented, elongate probe aperture In ac- 
cordance with the invention; 

45 Figure 8 shows a probe aperture in accordance with the invention having a cross-section of an elliptical 
shape of "width" 2 x f /, and "iengtti" 2 x f , ; 

Figure 9 is a plot in accordance with the invention of values of formation permeability versus preferred ra- 
tios of the radius of the impermeable pad to the radius of the probe aperture for laboratory testing with a 
mini-permeameter. 

so Figure 10 Is a table of values constructed In accordance with tihe Invention for an elliptic probe aperture 
having an aspect ratio of 0.2 oriented vertically and horizontally; 

Figure 11 is a graphic representation of the data presented in the first, second, third and sbcth columns of 
the table of Figure 10; 

Figure 12 is a graphic representation of ttie data presented in the first, fourth, fifth and sixtii columns of 
55 the table of Figure 1 0; 

Figure 13 is a table of values constructed in accordance with the invention for an elliptic probe aperture 
having an aspect ratio of 0.01 oriented vertically and horizontally; 

Figure 14 is a graphic representation of the data presented in the first, second, third and sixth columns of 
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the table of Figure 13; 

Figure 15 is a graphic represeutation of the data presented in the first fourth, fifth and sixth columns of 
the table of Figure 13; 

Figure 16 is a low chart of a preferred method for determining horizontal and/or vertical permeability in 
5 accordance with the invention; 

Figure 17 is a flow chart of a preferred method for determining horizontal and/or vertical permeability in 
accordance with the invention; 

Figure 18 is a table of values constructed in accordance with the invention for a rectangular probe aperture 
having an aspect ratio of 0.2 oriented vertically and horizontally; 
10 Figure 1 9 is a graphic representation of the data presented in the f irstz second third and sbcth columns of 
the table of Figure 18; 

Figure 20 is a graphic representation of the data presented in the first, fourth, fifth and sixth columns of 
the table of Figure 18; 

Figure 21 Is a table of values constructed In accordance with the Invention for a circular probe aperture 
IS and an elliptic probe aperture having an aspect ratio of 0.2 oriented horizontally; 

Figure 22 is a graphic representation of the data presented in the f irat. second, third and sixth columns of 
the table of Figure 21; and 

Figure 23 is a graphic representation of the data presented In the first, fourth, fifth and sbcth columns of 
the table of Figure 21. 

20 

DETAILED DESCRIPTION OF THE PREFERRED EI^BODIMENTS 

The Invention concerns nondestructive techniques for estimating the horizontal and/or vertical compo- 
nents of permeability of an anisotropic earth formation. As formations of interest typically comprise sedimentary 

25 rock, it is assumed that the formation Is isotropic In the horizontal directions, and has a smalier permeability 
in the vertical direction than in the horizontal. For purposes of this description, the "horizontal" directions are 
f hose generally parallel to the bedding planes of the rock, and the "vertical" direction Is generally perpendicular 
to the bedding planes of the rock. The term "formation" comprises a formation sample, such as a core plug 
taken from a borehole. In the case of a formation sample, "formation fluid" may be a liquid or a gas such as 

30 atmospheric air. It is noted that where a gas zone under consideration has been contaminated with liquid, meas- 
urements should be treated as if the formation sample is a liquid. 

In accordance with the invention, flow measurements are made to obtain values from which the perme- 
ability components of an earth formation are estimated. The flow measurements may be conducted in situ 
and/or in the laboratory using formation samples, in situ measurements are preferably made in a borehole with 

35 a formatk)n test tool having a probe aperture modified as described below. Formation test tools which may be 
employed include the Schlumberger RFT™ tester, MRTT^ tester and MDT™ tester. Laboratory measurements 
and measurements on outcrops are preferably made with a mini-permeameter having a probe aperture modi- 
fied as described below. 

The technique can be performed using a single probe. Pressure measurements are taken at the probe, 

40 through which fluid is forced to flow under subsfantially steady-state, single-phase conditions. For downhole 
measurements, the flow is preferably Induced by drawing formation fluid into the tool through the probe ("draw- 
down"). Alternately, fluid may be injected Into the formation through the probe ("injection"). Gas injection is pre- 
ferred for laboratory measurements with formation samples. Whether fluid is drawn into the probe or Injected 
out through the probe, a pressure disturbance is caused in the formation fluid. 

45 The technique may be used to determine permeability on a length scale sImOar to that of the Hassler core. 
Thus, permeability determined by this technique should be comparable to that obtained using the recognized 
standard procedure in the petroleum industry. 

Preferred methods of estimating horizontal and/or vertical penneabllity in accordance with the invention 
differ in at least two significant ways from the prior art methods described above. First, a probe having an aper- 

so ture of non-cirvutar cross-section is employed. The probe is that part of the tool or instrument In contact with 
the formation or formation specimen. Fluid is displaced through the probe aperture in making a measurement. 
The aperture is preferably shaped as a narrow slit, a smalt aspect ratio (width / length) being of more impor- 
tance than the exact shape of the cross-section. The slit shape allows f luhJ to be drawn or injected in a pattem 
which corresponds to the direction of measurement. For example, Figure 6 shows the probe oriented horizon- 

55 tally. As can be seen from the flow lines in Figure 6, the fluid enters the probe (in the case of draw>down) along 
the vertical axis V. Similarly. Figure 7 shows the ptobe oriented vertically. The flow lines in Figure 7 show the 
fluid entering the probe (in the case of a draw-down) along a horizontal axis X. The limit on the smallness of 
the aspect ratio results from a desire to avokJ clogging, and the size of the diameter (maximum length) of the 
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probe. The aspect ratio as defined (width/length) is less than 1.0. 

Second, measurenftents are taken during two pressure disturbances (e.g.. during two draw-downs), with 
the aperture oriented in two different directions with respect to the formation or formation specimen during the 
two measurements. For example, the aperture is oriented in a first direction (e.g., horizontal) during af irstdraw- 

5 down, and is oriented in a second direction (e.g., vertical) orthogonal to the first direction during a second draw- 
down. The "orientation" is the direction of the longest dimension of the aperture cross-section. 

A number of variations are possible. For example, the non-circular aperture cross-section may be generally 
elliptic or rectangular or of some other elongate or slit-like form. Instead of pressure draw-downs caused by 
withdrawal of fluid from the formation, pressure increases caused by injection of fluid into the formation may 

10 be used. A combination of a pressure draw-down and a pressure increase (injection) may be used in place of 
two draw-downs. Probes with two different aperture cross-sections may be used for the two pressure distur- 
bance (drawdown and/or injection) measurements - for example, one of the aperture cross-sections can be 
circular, provided the other aperture cross-section has a small aspect ratio (ratio of width to length). 

Determination of horizontal and/or vertical permeability in accordance with the preferred embodiments is 

15 based upon our derived relationship among the following parameters: the volumetric Dowrate, Q, and the vis- 
cosity, \i , of the fluid forced to pass through the aperture of the probe during draw-down or injection, the hor- 
izontal, kfj , and vertical, ky . components of the permeability of the formation, the pressure at the probe, Pp , 
the pressure of the formation far from the probe (equivalent to the pressure measured by the probe when the 
formation fluid is in its undisturt)ed state), P^ and the probe aperture dimensions 2 x f /, and 2 x f y. This rela- 

20 tionship Is obtained from the solution to the following boundary value problem wherein liquid is the fluid under 
constderatk)n: 



25 




P= Ppat|4 + Ir^landy = 0. (2) 

P-> P/as x2 + y2 + z^-^ 00 and y ^ 0 . (4) 
Due to the difference in compressibility between liquid and gas, the equations for gas become: 

35 
40 

= (Pp^atf^ + l^slandy = 0, (6) 

^ = 3^ = 0' (7) 

45 fP2;->(P^;asx2 + y2 + z^-^ooandygO. (8) 

P denotes the pressure field within the formation, and (x, y, z) denotes a rectangular Cartesian coordinate 
system oriented such that the x-axis and y-axis point in the horizontal directions and ttie z-axis points In the 
vertical direction, witii the y = 0 surface closely approximating the location of the borehole wall near the probe 
and the formation occupying the domain y ^ 0. In the case of the mini-permeameter, it is assumed that meas- 
^ urements are being made on a face of the formation sample which would satisfy tiiese conditions if it was still 
in tiie ground. For the moment, the cross-section of the probe aperture is assumed to have an elliptical shape 
of "widtti" 2 X Fft and length" 2 x f such as shown in Figure 8. (Examples of other possible aperture cross- 
sections are discussed below.) In Equation (1), the "V term, kJid^P/dz^, relates to formation permeability in 
ttie vertical directton and the %" term, ktid^P/dx^ ■»■ d^Pidy^, relates to formation permeability in an isotropic 
^ horizonfal plane Similarly, In Equation (5) the "V term relates to formation permeability in the vertical direction 
and the "V term relates to formation permeability in an isotropic horizontal plane. 
The desired relationship follows from the definition of volumetric flow rate, O, 
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(9) 



where Ap denotes the area of the aperture of the probe. The solution to this tx>undary-value problem appears 
in J.N. GCX)DIER et a/.. ELASTICITY AND PLASTICITY, John Wiley & Sons, inc., pp. 29-35 (1958). It gives 
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(10) 



where F denotes the complete elliptic integral of the first kind, and fp denotes the effective probe radius, defined 



20 as 



25 



30 



35 



40 



45 



SO 



= VV^TT^. (11) 



Ki4 and denote the dimensionless horizontal component and the dimensionless vertical component of the 
permeability, respectively. For liquid: 



For gas: 



(12) 



(13) 



It is assumed for in situ measurements that the aperture cross-section is sufficiently small compared to 
the radius of the hole (e.^., the well bore) containing the formation tester that the surface of the formation near 
the probe can be regarded as planar. For laboratory measurements (e.g., using a mini-permeameter and a for- 
mation sample), it is assumed that an impermeable pad surrounds ttie probe aperture to provide a hydraulic 
seal between the probe tip and the sample. The size of the pad and the formation sample are assumed to be 
large enough to justify the no-lux boundary condition on the entire y 0 surges (other than at the aperture) 
and the use of the semi-infinity domain (e.g., the "half-space" of D. GOGGIN et ai., A Thooretlcal and Expert 
mental Analysis of MInlpermeameter Response Including Gas Slippage and High Velocity Flow Effects, IN 
SITU, 12(1&2), pp. 79-116 (1988), at Figure 1). Figure 9 plots values of permeability, k, versus preferred ratios 
of Rpad/F^rDbe> where Rpad is the radius of the impermeable pad and Rp^be is the radius of the probe apertura. 
Pad dimensions for in situ measurement are less critical, in part due to the sealing effect of mud-cake at the 
borehole wall. 

The dimenstonless horizontal and vertical components of the permeability are determined as follows. Let 
2x2a and 2 x F/ denote the smallest and largest dimensions of the aperture of the probe, respectively. It will 
be recalled that we are Interested In any aperture having a small aspect ratio, i.e., the ratio ^Jii is a small num- 
ber. Avertlcal orientation of the probe aperture assumes f/, equals f^, and fy equals f/. A horizontal orientation 
of the probe assumes tfi equals f /, and f y equals f It is further assumed that two drawdowns are performed. 
During the first drawdown f iukJ flows through the probe at a volumetric f lowrate corresponding to Of, with the 
probe oriented vertically. During the second drawdown fluid lows through the probe at a volumetric flowrate 
corresponding to Q2, with the probe oriented horizontally. It is assumed that the values of Qf and are known; 
they need not be equal. This gives rise to the following two simultaneous equations containing only two un- 
knowns. Km and Kyp 
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The subscripts 1 and 2 refer to the pressure at the probe and flow rate through the probe corresponding 
to the first draw-down and the second draw-down, respectively, in the definitions of Kh and K^. The definition 
of the quantity Mfor liquid is given by: 

M = $^^l— g. (16) 

Ql(Pp2- Pf) 

The definition of the quantity Mfor gas Is given by: 



M 



(17) 



OW- Pf^Ppi' 

The value of quantity Mis readily obtained from the measured pressures and known flow rates, and is equivalent 
to both Khi/Kh2 and to Ky/f/Kyfz- The values of K^f and Kvt, hence the values of and are determined from 
the solution to the above set of equations. 

Values for K^i and K^t can be obtained by using a table such as Table 1 shown in Figure 10. The table is 
constructed from the above set of equations by evaluating the quantities M, K^^ K^i, Kh2 , and Kv2 over a range 
of values of the anisotropy, kt/ky, of the formation, for a given aperture aspect ratio . Table 1 is constructed 
foran elliptical aperture having aspect ratio i^i- 0.2 oriented vertically (subscript 1) and horizontally (subscript 
2). The equation makes use of the facts that kf/k^ = k^-^vi^ ^hz ^ ^hi^* and the value of f is known. That 
is, for a selected value of kt/ky^ equation ( 14) is used to evaluate Knf and equation (15) is used to evaluate 
Kh/M. The value of M is obtained by evaluating the ratio Khi/Kh2- Finally, Kvi and K^^are obtained by evaluating 
(ky/kn) X Khi and (kykt,) x Kh2. respectively. These evaluations determine a row In the table. Addittonal rows of 
the table are obtained by repeating these evaluations for the desired range of values for MCv^ 

To use the table, a value of M is calculated from measured pressure values and known flow rates of a set 
of pretest measurements made with the probe aperture oriented in the vertical direction during a first draw- 
down and in the horizontal direction during a second draw-down, or wee versa (see equation (16) for liquids 
and equation (17) for gases). The values of Khi and Kvi (or Kh2 and Ky2) in the same row as the calculated 
value of M represent the solution to the above set of equations. For example, if ijif equals 0.2 and M equals 
0.6732. then Table 1 (Figure 1 0) gives a value for Khi of 1 .905 and a value for K^ of 0.1 905. The exploit values 
of kt,, and k^ follow directly from the definitions of Khi and Kyi (or Kh2 and K^z), and the known values of Ppf, 

Ppi - (or Pp2. Pp2 - Pr). Qi (or Q2), n, and fp (/.a^ 

rp = + ). (18) 

Figures 11 and 12 graphically represent the data presented in Table 1 . In Figure 11 , the values of the ani- 
sotropy, k^h, and the dimensionless components of the permeability, Khi and Kvi are plotted versus values 
of calculated measurement factor M for an elliptic probe aperture having an aspect ratto of 0.2. The plotted 
values correspond to data presented in the first, second, third, and sixth colums of Table 1. The subscript 1 
denotes data characterizing the vertically oriented probe. In Figure 12, the values of the anisotropy. k^k^, and 
the dimensionless components of the permeability, Kh2 and Kv2, are plotted vereus values of calculated meas- 
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urement factor M for an elliptic probe aperture having an aspect ratio of 0.2. The plotted values correspond to 
data presented in the first, fourth, fifth, and sixth columns of Table 1. The subscript 2 denotes data character- 
izing the horizontally oriented probe. The values of the anisotropy. k^kn, and the dimensionless components 
of the permeabQity, K^i and Kyi (or Kh2 and Kyz), can be determined directly from these graphs. 

5 Table 2 (Figure 13) gives values for an elliptic aperture having an aspect ratio Fo/f/Of 0.01 oriented vertically 
and horizontally. The data of Table 2 is presented graphically In Figures 14 and 15. In Figure 14, the values of 
the anisotropy, k^n, and the dimensionless comronents of the permeability, Khi and Kyu are plotted veraus 
values of calculated measurement factor M for an elliptic probe aperture having an aspect ratio of 0.01 . The 
plotted values correspond to data presented in the first, second, third, and sixth columns of Table 2. The sub- 

10 script 1 denotes data characterizing the vertically oriented probe. In Figure 15, the values of the anisotropy, 
kyktf, and the dimensionless components of the permeability, Kh2 and Ky2* are plotted versus values of calcu- 
lated measurement factor M for an elliptic probe aperture having an aspect ratio of 0.01 . The plotted values 
correspond to data presented in the first, fourth, fifth, and sixth columns of Table 1. The subscript 2 denotes 
data characterizing the horizonfally oriented probe. The values of the anisotropy, k^^, and the dimensionless 

IS components of the permeability, Knt and Kyi (or and Kyz), can be determined directly from these graphs. 
It is also rather straight-forward to determine the propagation of error firom the measured quantity M to the 
predicfed quantifies k^ and ky. If there is a ±10% error in M, then the range of possible values of Kh and Ky 
corresponds to their values In rows bracl^eted by M equal to 1.1 x M and 0.9 x M, For example, if 2jti equals 
0.2 and M equals 0.67, then Table 1 (Figure 10) gives for the vertical probe 177 ^ Khi^ 2.07, or, 1.92 ± 7.6% 

20 error, and 0.97 ^ Kyi ^ 0.32, or, 0.21 ± 54% error, and for the horizontal probe 2.38 ^Kh2 = 3.42, or 2.90 ± 
17.8% error, and 0.17 ^ Ky^ ^ 0.44, or, 0.31 ± 43% error. In this case, the most accurate determination of Kh 
and Ky is obtained using the results from the vertical probe for Kh and the horizontal probe for Ky. 

When the aspect ratio of the probe aperture decreases in value, the error propagated also decreases. For 
examples, if tjti equals 0.01 and M equals 0.47 (corresponding to the same anisotropy as in the previous ex- 

25 ample), then Table 2 (Figure 1 3) gives for the vertical probe 3.1 5 ^ Khy ^ 3.30, or, 3.22 ± 2.2% error, and 0.24 
^ Kyi ^ 0.43, or. 0.33 ± 29.3% error, and for the horizontal probe 8.20 £ ^ 7.88, or 7.04 ±11.9% error, 
and 0.56 ^ Kv2 ^ 0.83, or, 0.70 ± 19% error. Again, the most accurate determination of Kh and Ky consists of 
using the results from the vertical probe for Kh and the horizontal probe for Ky, Note the improvemenf In accu- 
racy by using a probe with a smaller aspect ratio. 

30 Flow charts of preferred methods in accordance with the invention are given in Figures 16 and 17. The 
probe Is applied to the formation (or formation sample) with the aperture oriented In a first direction, preferably 
either horizontal or vertical (step 1610). The formation pressure is measured at the probe (step 1620). Muid is 
displaced through the probe for a first time period at a flow rate Qi (step 1630). Pressure at the probe Is meas- 
ured at the end of the first time period (step 1640). The probe is then withdrawn, rotated 90'', and reapplied to 

35 the formation (step 1650). Fluid Is displaced through the probe for a second time period at a flow rate (step 
1660). Pressure at the probe is measured at the end of the second time period (step 1670). Viscosity of the 
fluid is measured (step 1680). Values of horizontal permeability kh and/or k^, are determined from the aperture 
dimensions, the measured pressures, the low rates, and the fluid viscosity. 

A preferred embodiment of determining horizontal and/or vertical permeability values {e.g., of performing 

40 step 1690) is shown in Figure 17. Values are obtained for the aperture dimensions, the measured pressures, 
the flow rates, and the fluid viscosity, such as by the method of Figure 16 (step 1 710). A value for measurement 
factor M is calculated using the measured pressures and the flow rates (step 1720). Permeability factors Khi 
and Kyi (or Kh2 and Ky^) are evaluated using the aperture dimensions and the value for measurement factor 
M (step 1730). Values of k^ and/or ky, are determined from the permeability foctore, the aperture dimensions, 

45 the measured pressures, one or both of the flow rates, Qi and Q2, and the fluid viscosity. 

The steps of Figures 16 and 17 need not be carried out in the precise order given. For example, the for- 
mation pressure may be measured at the probe at any suitable stage in the process, or may be measured at 
aseparate probe. The viscosity of the displaced fluid may be determined at any time prior to determining values 
for kft and/or ky, by testing of a sample or by estimation or of herwise. 

so Other aperfure shapes may be used, such as that of a rectangle. For this case an approximate solution to 
the boundary value problem has been obtained. Instead of assuming that the pressure of the fluid takes on a 
constant value at the aperture, it is assumed that the velocity of the fluid leaving the formation is the same at 
every point of the aperture. Expressions have been derived relating 0, ^, kf^, ky, and Pp - Pf for the probe oriented 
both vertically and horizontally with respecf to the formation (formation sample) with an aperture having di- 

55 mansions 2 x f s by 2 x f /, where Pp denotes the average pressure over the aperture (see H.S. CARSI^W et 
ai., CONDUCTION OF HEAT IN SOLIDS, Oxford Science Publications (1959)). They are 
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where the definitions of Km, Kvi^ Khz, Kv2, Mand fp are the same as in the case gI the elllptlcaliy-shaped aper- 
ture, with the exception that Pp - Pf takes the piace of Pp-Pf. For liquids: 



(21) 
(22) 



and 



For gases: 



and 



Qi(Pa2- Pr) 



(23) 



(24) 
(25) 



Simuitaneous equations (19) and (20) can be solved using the same technique as before. For example, 
the variables M, K^u Kyu Khz, and have been evaluated over a range of values of ktJky for a rectangular 
aperture with aspect ratio equal to 0.2. The data are given in Table 3 of Figure 18, and presented graphically 
in Figures 19 and 20. Note the similarity between Table 1 (Figure 10) and Table 3 (Figure 18). 

Probe apertures of different shapes may be used for the two pressure disturbance measurements (e.g., 
draw-downs). One of the two probe apertures may be circular. For example, assume that probe 1 has a circular 
aperture of radius rpi and that probe 2 has an elliptical aperture of known aspect ratio F^/f / oriented horizontally 
with respect to the formation (or formation sample). The relevant relationships follow from the results for the 
elliptical aperture. They are 
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For gases: 
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. 4r^[(Pn2^- P/^P^/> . 

The value fp^, for the elliptical aperture is given by: 

r,2 = VF,2y2 + (33) 

A solution to simultaneous equations 27 and 28 can be obtained using the same method as described in 
the aboye examples. Table 4 of Figure 21 contains evaluations of M, K^, Kvt, K^z, and Ky? over a range of 
values of kt/ky for the case of a circular aperture and a horizontal elliptical aperture with aspect ratio tjtf equal 
to 0.2. These results are lllustFated graphically In Figures 22 and 23. 

While the foregoing describes and illustrates particular preferred embodiments of the invention, it will be 
understood that many modifications may be made without departing from ttie spirit of the invention. For exanv 
ple, It may be possible to use a first elongate shaped probe having width 2 x 2^1 and length Then, dig the 
second sampling In an orthogonal, second direction, a second elongate probe having width 2 x f .2 and length 
2 a is used. The two probes may differ in their overall dimensions. However, the mathematical Interpretafion is 
equivalent. The preferred embodiment presumes that the dimensions are the same for simplicity. Also, it may 
be possible to have a rectangular shaped probe instead of the elliptical shaped probe during tiie second sam- 
pling, whBe having a circular probe during the first sampling or Wee versa. We Intend ttie following claims to 
cover any such modifcations as fell within the true spirit and scope of the invention. 



Claims 

40 

1. A method of estimating permeability of an et formation in at least one of two orthogonal directions, the 
formation containing a formation fluid, comprising tiie steps of: 

a. measuring a pressure P^ of the formation fluid; 

b. creating a pressure disturbance in the formation fluid by displacing fuid through a probe aperture for 
46 a first time period at a first flow rate Of, tiie probe aperture having an elongate cross-section of width 

2 X f « and length 2 x F/ and being oriented in a first direction; 

c. measuring a pressure Ppf of the fluid substantially at the end of tiie first time period; 

d. creating a pressure disturbance in the formation fluid by displacing fluid through a probe aperture 
for a second The period at a second rate tiie probe aperture having an elongate cross-section of 

GO widtti 2 x F« and lengtti 2 x f / and being oriented In a second direction orttiogonai to said first direction; 

e. measuring a pressure Pp2 of the fluid substantially at the end of tiie second time period; 

f. determining a value \i for viscosity of fluid in the formation; and 

g. determining a value of permeability in at least one of said first and second directions from the aperture 
widtti 2x f«and ttieaperturelengtti 2x , the measured pressure Pf, at least one of measured pressures 

55 Ppf and Pp2, at least one of ttie flow rates Of and Q^. and the determined value \i for viscosity of fluid 

in the formation. 

2. The mettiod of daim 1. wherein step g. comprises ttie steps of: 
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i. calculating a measurement factor M from the measured pressures Ppi and Pp2 and form the flow 
rate Of and Q2 ; 

II. determining a value of dimenslonless quantity K)^/ representative of the horizontal permeability of the 
formation and a value of a dimenslonless quantity representative of the vertical permeability of the 
formation, based on the calculated measurement factor M and the aperture length 2 x f « and the aper- 
ture lenght 2 x f j ; and 

Hi. determining a horizontal permeability value from the values of quantity Khh the aperture width 2 
X and the aperture length 2 x the measured pressure P/, at least one of measured pressures Ppi 
and Pp2 , at least one of flow rates Qi and O2, and the determined value \i for viscosity of fluid in the 
formation. 



3. The method of daim 1, wherein step g. comprises the steps of: 

I. calculating a measurement factor M from the measured pressures P^ Ppi and Pp2 and from the flow 
rates Qi and Q2 in accordance with the relationship 
IS w - Q2(P/>i - Pd . 

Qi(Pp2-Pf)' 

il. determining a value of a dimenslonless quantity Khi representative of the horizontal permeability of 
the formation and a value of a dimenslonless quantity K^i representative of the vertical permeability of 
the formation, based on the calculated measurement factor M and the aperture width 2 x f ^ and the 
20 aperture length 2 x f / In accordance with the relationships 
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where F denotes the complete elliptic integral of the first kind; 

Hi. determining a horizontal permeability value kh from the values of a quantity Kh/ comprising one of 
the quantities Khi and K^f/M. the aperture width 2 x f « and the aperture length 2 x F/, the measured 
pressure Pf, a measured pressure Pf^ comprising one of measured pressures Ppi and Pp^, a flow rate 
Q„ comprising one of flow rates Qf and Q2, and the determined value ^ for viscosity of fluid in the for- 
mation in accordance with the relationships 

ix _ 4rJiPax - Pdk„ 
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The method of daim 1, wherein step g. comprises the steps of: 

L calculating a measurement factor M from the measured pressures Ppf and Pp2 and from the flow 
rates Q| and Q2 in accordance with the relationship 

w . Q2(Poi- Ph . 
Qi(Pp2-Pr)' 

il. determining a value of a dimensionless quantity Khj representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity Kyn representative of the vertical permeability 
of the formation based on the calculated measurement factor M and the aperture width 2 x F« and the 
aperture length 2 x f/ in accordance with the relationships 
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where F denotes the complete elliptic integral of the first kind; 

iii. determining a horizontal permeability value kf, from the values of quantity Khu the aperture with 2 x 
f s and the aperture length 2 x f /, the measured pressure Pf, the measured pressure Ppt, the flow rate 
Qi, and the determined value ^ for viscosity of fluid in the formation in accordance with the relationships 

^ 4rp(Pjt - Pf)kh 



^H1 



fp = ye^^zVJ^ , 

The method of daim 1, wherein step g. comprises the steps of: 
i. calculating a measurement factor M from the measured pressures Pf, Ppi and Pp2 and from the flow 
rates Of and O2 in accordance with the relationship 

j^, Q2(Poi- Pf) . 
Q^(Pp2'P^' 

11. determining a value of a dimensionless quantity Khi representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity K^f representative of the vertical permeability 
of the formation, based on the calculated measurement factor M and tiie aperture widtti 2 x f « and the 
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aperture length 2 x in accordance with the relationships 
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where F denotes the complete elliptic integral of the first kind; 

iii. determining a vertical permeability vaiue from the values of a quantity comprising one of quan> 
titles Kyi and Kv/M» the apeture width 2 x f « and the aperture length 2 x f /, the measured pressure Pf, 
a measured pressure Pp; comprising one of measured pressures Ppt and Pp2, a flow rate Q„ comprising 
one of flow rates Qf and Qz, and the determined value ^ for viscosity of fluid in the formation in accor- 
dance with the relationships 

,2>2 + if /I . 



6. The method of dalm 1. wherein step g. comprises the steps of: 

i. calculating a measurement factor Af from the measured pressures Pf, Ppi and Pp; and from the flow 
rates Q| and in accordance with the relationship 

QyiP^^- Pf)Ppt' 

ii. determining a value of a dimensionless quantity representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity Ky^ representative of the vertical permeability 
of the formation, based on the calculated measurement factor M and the aperture width 2 x f « and the 
aperture length 2 x f / In accordance with the relationships 
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where F denotes the complete elliptic integral of the first kind; 

ill. determining a horizontal permeability value from the values of a quantity Kh/ comprising one of 
the quantities Khi and KHf/M, the aperture width 2 x f « and the aperture length 2 x f/, the measured 
pressure Pf, a measured pressure Pf^ comprising one of measured pressures Pp, and Pp2, a flow rate 
Qn comprising one of flow rates Qi and Q2. and the determined value ^ for viscosity of fluid in the for- 
mation in accordance with the relationships 

^ . 4rJ(P,j^- P^pjk, 

Tp = ,2>2 + ef/2 . 

The method of daim 1, wherein step g. comprises the steps of: 

i. calculating a measurement factor M from the measured pressures Pf, Ppi and Pp2 and from the flow 
rates Of and Q2 in accordance with the relationship 

ii. determining a value of a dimensionless quantity Km representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity Ky^ representative of the vertical permeability 
of the formation, based on the calculated measurement factor M and the aperture width 2 x f « and the 
aperture length 2 x f / in accordance with the relationships 
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where F denotes the complete elliptic integral of the first kind; 

iil. determining a vertical permeability value ky from the values of quantity Kyr comprising one of quan- 
tities Kvi and K^i/M, the aperture width 2xig and the aperture length 2 x f/, the measured pressure 
Pf, a measured pressure Pf^ comprising one of measured pressures P^i and Pp2, a flow rate Q„ com- 
prising one of flow rates Qi and Q2, and the determined value \i for viscosity of fluid in the formation 
in accordance with the relationships 

. 4rJ(P^/2- P?y2PplK 



8. A method of estimating permeability of an earth formation in at least one of two orttiogonal directions, 
comprising the steps of: 

a. measuring a pressure of fluid in the formation; 

b. creating a pressure disturbance in the formation fluid by displacing fluid through a probe aperture 
for a first time period at a first flow rate Qf , tiie probe aperture having an elongate cross-section of width 
2 X ^0 and length 2 x f / and being oriented in a first direction; _ 

c. measuring pressure of ttie fluid substantially at ttie end of the first period to obtain a value Ppi of 
average pressure over the aperture; 

d. creating a pressure disturbance in the formation fluid by displacing fluid through a probe aperture 
for a second time period at a second rate Q2, the probe aperture having an elongate cross-section of 
widtii 2 X f « and length 2 x and being oriented in a second direction orthogonal to said first direction; 

e. measuring pressure of ttie fluid substantially at the end of the second time period to obtain a value 
Pp2 of average pressure over ttie aperture; 

f. determining a value \i for viscosity of fluid in the formation; and 

g. determining a value of permeability in at least one of two orthogonal directions from the aperture 
widtii 2 X f 5 and ttie aperture iengtti 2 x the measured pressure Pf, at least one of the average pressure 
values Ppf and Pp2. at least one of ttie flow rates Q| and Q2, and ttie determined value \i for viscosity 
of fluid in the formation. 
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9. The method of daim 8, wherein step g. comprises the steps of: 

l^calculab'ng a measurememt factor M from the measured pressure Pf, the average pressure values 
Ppt and Pp2. and the flow rates Qi and Q^; 

II. determining a value of a dimenslonless quantity representative of the horizontal permeabOlty of 
5 the formation and a value of a dimenslonless quantity Ky representative of the vertical permeability of 

the formation, based on the calculated measurement factor M and the aperture dimensions 2 x f « and 
2xf,; 

ill. determining a horizontal permeability value kf, from the values of quantity Knh the aperture^imen- 
slons 2 X f 8 and 2 x f /, the measured pressure Pf. at least one of the average pressure values Ppf and 
10 Pp2p at least one of the flow rates Qf and Q^, and the determined value ^ for viscosity of fluid In the 

formation. 



IS 



20 



10. The method of dalm 8, wherein step g. comprises the steps of: 

i. calculating a measurement factor M from the measured pressure P^ the average pressure values 
Ppf and Pp2i and the flow rates and Q^; 

ii. determining a value of a dimenslonless quantity Kh representative of the horizontal permeability of 
the formation and a value of a dimenslonless quantity representative of the vertical permeability of 
the formation, based on the calculated measurement factor M and the aperture dimensions 2 x f , and 

2xf;; 

ill. determining a vertical permeability value /c^ from the values of quantity Kyr, the aperture^imensions 
2 x f « and 2 x f ^ the measured pressure P/v at least one of the average pressure values Ppf and Pp2, 
at least one of the flow rates Qf and Q2, and the determined value \i for viscosity of fluid in the formation. 
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11. The method of dalm 8, wherein step g. comprises the steps of: ^ 
i. calculating a measurement factor Mfrom the measured pressure Pf, the average pressure values Ppf 
and Pp2, and the flow rates Qf and Q2 in accordance with the relationship 

^ . Q2(P|l1- Pf) . 

Qi(Pp2- Pd' 

11. determining a value of a dimenslonless quantity Km representative of the horizontal permeability of 
the formation and a value of a dimenslonless quantity Kyi representative of the vertical permeability 
of the formation, based on the calculated measurement factor M and the aperture dimensions 2 x f « 
and 2 X f / in accordance with the relationships 
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iii. determining a horizontal permeability value from the value of a quantity comprising one of the 
values Khi and Knf/M, the aperture dimensions 2 x f « and2 x f the measured pressure Pf, a pressure 
value comprising one of the average pressure values Ppi and a flow rate Q„ comprising one of 
the flow rates Qi and Q^, and the determined value \x for viscosity of fluid in the formation in accordance 
with the relationships 
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12. The method of daim 8, wherein step g. comprises the steps of: 

i. calculating a measurement factor Mfrom the measured pressure the average pressure values Ppy 
and Pp2. and the flow rates Qi and in accordance with the relationship 

ii. determining a value of a dimensionless quantity K^^ representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity Ky^ representative of the vertical permeability 
of the formation, based on the calculated measurement factor M and the aperture dimensions 2 x 
and 2 X f / in accordance with the relationships 
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13. 



III. determining a vertical permeability value /r^ from the values of a quantity Kv7 comprising one of quan- 
tities Ky/i and Ky/i/M, the aperture dimensions 2 x and 2 x thj measured pressure P^ a pressure 
value P^ comprising one of the average pressure values Ppf and Pp2, a flow rate Q„ comprising one of 
the flow rates Qf and O2. and the determined value ^ for viscosity of fluid in the formation in accordance 
with the relationships 

I, _ 4rp(Pp/ - P^y 

The method of daim 8, wherein step g. comprises the steps of: 

i. calculating a measurement factor Mfrom the measured pressure Pf, the average pressure values Pp^ 
and Pp2, and the flow rates Of and Q2 in accordance with the relationship 
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ii. determining a value of a dimensionless quantity Khi representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity Kyi representative of the vertical permeability 
of the formation, based on the calculated measurement factor M and the aperture dimensions 2 x f< 
and 2 X f f In accordance with the relationships 
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ili. determining a horizontal permeability value from the value of a quantity Khi comprising one of val- 
ues Khi and Khi/M, the aperture dimensions 2 x f« and 2 x the measured pressure a pressure 
value Ppi comprising one of the average pressure values Ppi and Pp^, a flow rate comprising one of 
the flow rates Qi and Q2, and the determined value \i for viscosity of fluid In the formation in accordance 
with the relationships 
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14. The method of daim 8, wherein step g. comprises the steps of: _ 

i. calojlatlng a measurement factor Mfirom the measured pressure Pf, the average pressure values Ppi 
and Pp2, and the flow rates Qi and Qz in accordance with relationship 

ii. determining a value of a dimensionless quantity Km representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity Ky^ representative of the vertical permeability 
of the formation, based on the calculated measurement factor M and the aperture dimensions 2 x f « 
and 2 X f/ in accordance with the relationships 
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iii. determining a vertical permeability value k^ftom the value of a quantity Ky/t comprising one of values 
Kvi and Ky^/M, the aperture dimensions 2 x and 2x f,, thejneasured pressure P^, a pressure value 
comprising one of the average pressure values Ppf and Pp2, a flow rata comprising one of the 
flow rates Qf and Q2, and the determined value ^ for viscosity of fluid in the formation in accordance 
with the relationships 

4rJ(Pp/- P^pJ^v 
= + . 



Kv7 



45 



15. A method of estimating permeability of an earth formation in at least one of the horizontal and vertical 
directions, the formation containing a formation fluid, comprising the steps of: 

a. measuring a pressure Pf of the formation fluid; 

b. creating a pressure disturbance in the formation fluid by displacing fluid through a first probe aperture 
for a first time period at a first fow rate Of, the first probe aperture having a circular cross-section of 
radius 

c. measuring a pressure Ppi of the fluid substantially at the end of tiie first time period; 

d. creating a pressure disturbance in the formation fluid by displacing fluid through a second probe aper- 
ture for a second time period at a second rate Q2, the second probe aperture having an elongate cross 
section of width 2 x f « and length 2 x f /; 

e. measuring a pressure Pp2 of the fluid substantially at the end of the second time period; 

f. determining a value ^ for viscosity of fluid in the fornrtation; and 

g. determining a value of permeability In at least one of the horizontal and vertical directions from the 
aperture dimensions 2 x f 2 x f / and r^^ the measured pressure Pu at least one of the measured pres- 
sures Ppf and Pp2, at least one of the flow rates and Q2, and the determined value \k for viscosity of 
fluid in ttie formation. 



16. The method of daim 15, wherein step g. comprises the steps of: 

i. calculating a measurement factor M from the measured pressures Pf, Ppi and Pp; and from the flow 
^ rates Qf and Q^; 

11. determining a value of a dimenslonless quantity K^/ representative of the horizontal permeability of 
the formation and a value of a dimenslonless quantity representative of the vertical permeability of 
the formation, based on the calculated measurement factor Aland the aperture dimensions 2 x f^. 2 x 
and Tpf; and 

^ iii. determining a horizontal permeability value /t/, from the values of quantity an aperture dimension 

fpm comprising one of values fpf and rp2 where rp2 Is a function of 2 x f « and 2 x f 1; ttie measured pressure 
P/; at least one of measured pressures Ppf and Pp2; at least one of flow rates Q| and Q2; and tiie de- 
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termined value \i for viscosity of fluid In the formation. 



17. The niethod of daim 15.whereln step g. comprises the steps of: 

i. calculating a measurement factor M from the measured pressures Ppi and Pp2 and from the flow 
5 rates Of and 

11. determining a value of a dimensionless quantity Kh representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity Kw representative of the vertical permeability of 
the formation, based on the calculated measurement factor Aland the aperture dimensions 2 x 2 x 
f/andTpf! and 

10 iii- determining a vertical permeability value ky from the values of quantity Kw ; an aperture dimension 

fpm comprising one of values r^i and where rp2 is a function of 2 x f « and 2 x the measured pressure 
P/i at least one of measured pressures Ppi and Pps; at least one of flow rates Qi and Q2; and the de- 
termined value \i for viscosity of fluid In the formation. 

18. The method of daim 15, wherein step g. comprises the steps of: 
i. calculating a measurement factor M from the measured pressures P^ Ppi and Pp2 and from the flow 
rates Q| and O2 in accordance with the relationship 

Q2(P.i- Pf) . 

Ql(Pp2- Pf)' 

11. determining a value of a dimensionless quantity Khi representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity Kyi representative of the vertical permeability 
of the formation, based on the calculated measurement factor M and the aperture dimensions 2 x f s* 
2 X f / and fpi in accordance with the relationships 
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where F denotes the complete elliptic Integral of the first kind; 

ill. determining a horizontal permeability value ^/,from the value of a quantity comprising one of quan- 
tities Khi and Kh/M. a value Tp^ comprising one of values fpi and rp2, the measured pressure Pf, a meas- 
50 ured pressure P^ comprising one of measured pressures Ppi and Pp2, a flow rate Qn comprising one 

of flow rates Qf and Q2, and the determined value ^ for viscosity of fluid in the formation In accordance 
with the relationship 

^ , 4W(P^- P,yc, 

55 

19. The method of daim 15, wherein step g. comprises the steps of: 

1. calculating a measurement factor Af from the measured pressures P^ Ppi and Pp2 and from the flow 
rates Q| and Q2 In accordance with the relationship 
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ii. determining a value of a dimensionless quantity representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity Kvt representative of the vertical permeability 
of the formation, based on the calculated measurement factor M and the aperture dimensions 2 x f«, 
2 X f f and Tpf In accordance with the relationships 
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where F denotes the complete elliptic Integral of the first kind; 

ill. determining a vertical permeability value /r^ from the values of a quantity Kvr comprising one of the 
values Kvi and Kvf/M, a value fpm comprising one of values fp, and rp2. the measured pressure Pf, a 
measured pressure Pp/ comprising one of measured pressures Ppi and Pp2, a flow rate 0^ comprising 
one of flow rates Qi and Q^ and the determined value \i for viscosity of fluid In the formation in accor- 
dance with the relationship 

20. The method of dalm 15, wherein step g. comprises the steps of: 

I. calculating a measurement factor M from the measured pressures Pf, Ppi and Pp2 and from the flow 
rates Qi and in accordance with the relationship 

Qi(V- Wi' 

II. determining a value of a dimensionless quantity Khi representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity Kvi representative of the vertical permeability 
of the formation, based on the calculated measurement factor M and the aperture dimensions f «, f / and 
Tpt in accordance with the relationships 
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where F denotes the complete elliptic Integral of the first kind; 

ill. determining a horizontal permeability value ki, from the values of a quantity Knt comprising one of 
quantities Khi and Khi/M, a value r^m comprising one of values rpi and Tp^, the measured pressure Pf, 
a measured pressure Pp; comprising one of measured pressures Ppf and Pp2* a flow rate Q„ comprising 
one of flow rates Qt and Q^, and the determined value \i for viscosity of fluid In the formation In accor- 
dance with the relationship 

pressure P/^ the measured pressure Ppt, the flow rate Qu and the determined value p for viscosity of 
fluid In the formation In accordance with the relationship 

4rpi[(P^2- P^Ppil/g, 
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21. The method of claim 15, wherein step g. comprises the steps of: 

i. calculating a measurement factor Mfrom the measured pressures Pf, Ppi and Pp2 and from the flow 
rate Qi and Qz In accordance with the relationship 

ii. determining a value of a dimensionless quantity Knt representative of the horizontal permeability of 
the formation and a value of a dimensionless quantity K^i representative of the vertical permeability 
of the formation, based on tiie calculated measurement factor M and the aperture dimensions 2 x f „ 
2 X ^/ and Tpt in accordance with the relationships 
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where F denotes the complete elliptic integral of the first kind; 

ill. determining a vertical permeability value ky from the values of a quantity Kyi comprising one of the 
values Kyi and Kyi/M, a value fp^ comprising one of values r^i and rp2, the measured pressure Pf, a 
measured pressure Ppj comprising one of measured pressures Ppf and Pp2, a flow rate Q„ comprising 
one of flow rates and Qj, and the determined value ^ for viscosity of fluid in the formation In accor- 
dance with the relationship 
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